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Abstract 

We revisit the celebrated Peierls-Onsager substitution employing the magnetic pseudo- 
differential calculus for weak magnetic fields with no spatial decay conditions, when the non¬ 
magnetic symbols have a certain spatial periodicity. We show in great generality that the 
symbol of the magnetic band Hamiltonian admits a convergent expansion. Moreover, if the 
non-magnetic band Hamiltonian admits a localized composite Wannier basis, we show that 
the magnetic band Hamiltonian is unitarily equivalent to a Hofstadter-like magnetic matrix. 
In addition, if the magnetic field perturbation is slowly variable, then the spectrum of this 
matrix is close to the spectrum of a Weyl quantized, minimally coupled symbol. 
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1 Introduction and main results 


The mathematical description of the quantum theory of solids is mainly based on the spectral and 
dynamical analysis of Schrddinger-type operators with periodic coefficients. A central problem is 
concerned with the response of a periodic system when submitted to an external magnetic field 
which does not vanish at infinity. The main technical difficulty in dealing with such a problem 
comes from the fact that the perturbation represented by the magnetic field is no longer an analytic 
perturbation. In particular, the difference between the perturbed and the unperturbed observables 
is not small in operator norm for small values of the magnetic field intensity and a norm convergent 
perturbative approach is impossible. 

The solution proposed by physicists (see section 7.3 in [51]), the so-called Peierls-Onsager 
substitution, is to isolate as main contribution the magnetic quantization of the Bloch eigenvalues 
and, modulo some gauge transformation, to obtain a power series development with respect to 
the magnetic field intensity. In that framework, magnetic quantization means just to replace the 
canonical moment variables with the magnetic ones by the so called minimal coupling procedure. In 
mathematical terms this amounts to give a precise meaning to functions of the operator 0 —A(iVe) 
for 9 G Tf the d-dimensional torus in the momentum space (see m and references therein). 
Moreover, one has to to deal with an apparent lack of gauge covariance of this procedure. Although 
a rather rich literature is devoted to the precise mathematical treatment of this problem (see 
m n 1 nn na da m nil im ) , a complete understanding is still missing. 

The main objective of our paper is to show that the magnetic pseudodifferential calculus de¬ 
veloped in na dzi dH] and its integral kernel counterpart developed in |30j can give a precise 
mathematical framework and a gauge covariant formulation of the above physical proposal by 
-roughly speaking- replacing Effk — A{x)) of the Peierls-Onsager substitution method (here En is 
the n-th Bloch eigenvalue of the unperturbed system) with the magnetic pseudodifferential oper¬ 
ator Op^{En). In this paper we shall mainly concentrate on the so-called isolated spectral island 
situation (see Hypothesis 11.51) and leave the case without a spectral gap for a forthcoming paper. 

Let us point out from the very beginning that our treatment is gauge covariant. The magnetic 
fields are supposed to be smooth and bounded but no condition of slow variation is imposed in 
general. Of course, under an extra hypothesis of slow variation, stronger results may be obtained 
and we shall compare them with those already existing in the literature. Let us also underline that 
the results obtained cover a large class of Hamiltonians described by periodic pseudodifferential 
operators. 

One of the conclusions of our analysis is that once we can isolate a spectral island of a periodic 
Hamiltonian and define symbols of the associated band operators (Hamiltonian, projection, etc.), 
then their corresponding magnetic counterparts, while being singular perturbations of the free 
ones, can be well approximated in the norm topology by the magnetic quantization of the ‘free’ 
symbols. Moreover, when a composite Wannier basis is supposed to exist for the unperturbed 
spectral projection mm, then a generalization of the results in m is obtained. 

1.1 The framework and some notation 

We shall consider a d-dimensional configuration space with d > 2 and use the notation A = 
Although its dual is canonically isomorphic to we shall prefer to use the notation X* in 
order to emphasize the dual variable. We shall denote by (•,•) : X* x A —)> R the duality 
relation. We shall consider S := A x A* as a symplectic space with the canonical symplectic 
form a°{X,Y) := {ffy) — {rj,x) where X = and Y = {y,r]). 

We shall consider an algebraic basis {ei,..., e^} of R^^ and the lattice T = it gives 

an injective homomorphism lA G X that allows us to view T as a discrete subgroup of A. We 
consider the quotient group A/T that is canonically isomorphic to the d-dimensional torus T'^ that 
we shall denote by T when the subgroup T is evident; let us denote by tt : A —T the canonical 


2 


projection onto the quotient. Let us fix an elementary cell: 


E := 


y = 


d 

i=i 


-1/2 <t, <1/2, Vje{l,...,d} 


and the unique decomposition it induces for any x € X: x = [x\+ x with [a;] G T and x € E. The 
dual lattice of T is defined as 


r, := {7* e T”* I ( 7 *, 7 )/( 27 r) G Z, V7 S T} . 

Considering the dual basis {e^,..., e^} C X* defined by (e*, Cfc) = {2TT)Sjk, we have T* = 0^^;^Ze*. 
We define T* := <Y*/r* and S* by a definition similar to that of E but with respect to the dual 
basis {ej}i<j<d- We notice that T* is isomorphic to the dual group of T and T* is isomorphic to 
the dual group of T (in the sense of abelian locally compact groups); moreover let tt* : X* —>• T* 
be the canonical quotient projection. For any 7 G F, 7* G F* or a; G fb, and ^ G X* we shall 
denote by Ux and uj the corresponding characters on T*, T, lb and lb* respectively, i.e. 

a^{B) := a^.(^(x)) := a,(^) = a^{x) := 

Given some finite dimensional real Euclidean vector space V, we shall denote by ^(V) and 
^'(V) the Schwartz space of test functions and resp. its dual, the space of tempered distributions, 
on V; we shall denote by (•, •)v : o5^'(V) x J^(V) —>■ C the canonical duality. We shall constantly 
use the notation < v >:= \/\ + |up for any u G V. We shall consider the spaces BCiV) of 
bounded continuous functions with the || • ||oo norm. We shall denote by C°°{y) the space of 
smooth functions on V and by G^/V) and by BC°°{y) its subspace of smooth functions that are 
polynomially bounded together with all their derivatives or smooth and bounded together with 
all their derivatives. 

Let us recall a class of Hormander type symbols on S that we shall use. For any s G K and 
any p G [0,1] we denote by 

S;{X) := {E G G“(S) I < (X),V(a,&) G x N'"} ( 1 . 1 ) 


where v"S{f) := 


sup , V(a,6) G X N'^; we recall that |6| := f) bj. 

(x,{)eH j=i 

A symbol F in Sp{X) is called elliptic if there exist two positive constants R and C such that 
|G(a:,^)| > G(^)® for any (x,^) G S with |^| > R. 

Let us recall that for h G 5'™(A’) the Weyl quantization associates the operator Op{h) defined 
on S^{X) (see formula (13.31) with A = 0), and having a natural extension by duality to J^'(fb). 

We denote by t the action by translations both of X and of X* on the tempered distributions 
defined on the respective spaces, i.e. [tzF){x) := F{x + z) for any F G C{X), and any {z,x) G 
X X X. We shall consider symbols of Hormander type SY^{X) satisfying ((r.y ( 8 > i)h) {x, ^) = h{x, 
for any 7 G F and for any (x,^) G S and we denote by 5'™(A’)r the space of these symbols. We 
shall use the following notations: 


• y^{X) := {u G y'{X) I TjU = u, V7 G F}, the space of F-periodic distributions on X. 

• y{T) := G°°(T) with the usual Frechet topology; y'{T) is the dual of y{T). We shall 
denote by (•, •)t the natural bilinear map defined by the duality relation on y'(T) x y{T). 

Remark 1.1. We have the identification y{T) = y^{X) fl C°°{X) obtained by transporting a 
F-periodic function to the quotient by F. It is also well known that the spaces .5^/(A’) and y’(T) 
can be identified by a natural topological isomorphism. 

Hypothesis 1.2. From now on we shall consider a fixed real elliptic symbol h G 5'™(A’)r for 
some m > 0. 

Let us recall that for a h G S'^{X) which is real, elliptic, and with m > 0, the operator Op{h) 
has a self-adjoint extension El with domain Jtf'^{X) := {/ G L^{X) \ (1 — A)"*/^/ G L^(A’)}. Let 
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US consider its resolvent {H — 3 )“^ at 3 G C \ cr{H); it is known that it is a pseudodifferential 
operator with a symbol r^{h) of class Let us notice that if h £ 5'™(<Y)r, the operator 

Dp(/i) commutes with the action of L through translations and thus we can define its restriction 
to test functions on the torus: 


( 1 . 2 ) 

This allows us to consider its self-adjoint extension Hf defined on fl 

and acting in the Hilbert space L^(T) = S^^{X) n Lf^^{X). 

Given a Hilbert space K, and two vectors u and v in /C, we shall denote by |u){u| the rank 1 
operator in K associated to it, i.e. the operator given by |u)(u| w := u(v,w)jc, for any w € 1C. 

We shall need the following notation for different types of Fourier transforms: 

Ts^(X) := {21,)-'^ ^ dY, G ^(S); := (^, x); 


Txm := ( 2 ^)-"/^ [ 
Jx 

TrviO) ■■= 

7Gr 


\f(t)es^{xy, 

VvGi\r). 


For any function T* 9 0 !->■ (f{9) G C we denote by (p{S,) := its periodic extension to X*. 

Finally let us recall the Poisson formula ([H]) that we shall use several times: 


\E.\ 


7er 


y~! ^7*> 

7 * 6 r, 


(1.3) 


as distributions in dX^X) with series converging in the weak sense. 


1.2 The Bloch-Floquet transform 

Let us recall some basic facts concerning the Bloch-Floquet method and fix some notation that 
will be used in the following. We consider 


^ {p G Ll^{X X X*) I t^F = a-fF V 7 G F,T.y.F = F V 7 * G F*} 


(1.4) 


with the Hilbertian norm 


F 




d^ dx and the the Bloch-Floquet unitary map 




{^Ti){x,i) = -j) 

7Gr 


with its inverse having the explicit form 


(%-^f) (xo+7 ) 


|T,r^ f a-,{9)Fixo,9)d9. 


(1.5) 


( 1 . 6 ) 


Remark 1.3. For any 7 G F we notice that = cr-y'^r- 

Let us briefly recall the direct integral structure of that is very useful in the study of periodic 
'LDO. Using Fubini Theorem we consider each function in as a function defined on T* and taking 
values in Lf^^{X) and define for each 0 G T* the complex linear space 

:= {/ € I = ay(d)/ = a,(7)/} . (1.7) 

For any 0 G T* we define the Hilbertian norm ||/||^j, = J^\f{x)ydx and we can define a family 
of unitary operators Yg : ^ L^(T) by the formula YbF{e{x)) := a-e{x)F{x,9). Following 
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the procedure in [S] one can consider on the field of Hilbert spaces over the compact 

smooth manifold T» the measurable structure defined by the constant field L^(T) and the field of 
unitaries {'fejeej, and define the associated direct integral of Hilbert spaces fj ^gdO that will 
be unitarily isomorphic to ^. 

For any / € ^{X) we can write 


(^rDp(h)/)(0,cr) = 5]e-*<^-'^>(Dp(h)/)(a;-7) 

7er 


7Gr 


(x), 


with the series converging for the weak topology of tempered distributions. Thus 


«p(h)'^r-i = r Dp{h)\^ dO. (1.8) 

7t* 


Definition 1.4. We consider the following spaces of tempered distributions on 

1. ^/(S) := {F G y'{E) I t^F = ct^F,^-/ G F, t^F = G F,} , 

2. := {F G ^/(S) I ((1 - A)*/2 (g) ]l)i^ G ^} , Vs G R 

3 . ■■= {fe^0\{t- A)«/ 2 / G , Vs G R, V 0 G T*. 


The operator H := acting in the Hilbert space ^ is self adjoint on the domain 

We notice that the operator takes any smooth function T* 9 0 !->■ f{6) G ^(m),e 

into the smooth function T, 9 0 n- Dp(/i)|^^ f[0) G ^g. We can also consider H{9) as the 
self-adjoint operator defined by Dp(/i)|^^ in the Hilbert space J^g with domain and due 

to its ellipticity we have that ^(rn),e = {/ S ^g \ (Dp(/i)|^^) / G ,^g}. On these subspaces the 


quadratic graph-norms 


I {m),g 




LHE) 


define a field of Hilbert spaces that are 


unitarily transformed into 'H'"(T) through the unitaries {Ve}egT, ■ Using once again the theory in 
[S] we can define the direct integral Hilbert space d9 that will be unitarily isomorphic to 

^(^rn) ■ Due to the above arguments, we can also define the measurable field of bounded operators 


^H{6) G and the associated direct integral 


H = ^rH%-^= / H{e)de G 
Jt* 

One has the identity a-g Op{h)\^^ ag = Dp{Tgh)\^^ as operators acting in L^(T). For the 
convenience of the reader we remind the following fundamental result: 

Bloch-Floquet Theorem. Under Hypothesis 11.21 the self-adjoint operator F[ has the following 
spectral properties: 

1. There exist a family of continuous functions T, 9 0 i-9 Xj{0) G R indexed by j G N* 
such that \j{6) < Xj+i{6) for every j G N* and 0 G T*, and 

a{Hie)) = y {A,(0)}. 

iGN* 


2. There exists a family of measurable functions T* 9 0 g->. G indexed by 

j G N* such that ||7(^)ll.^s = 1 

H{9)(t)j{9) = Xj{9)(l)j{9), Vj G N*, V6» G T,. 
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Hypothesis 1.5. There exists a bounded open interval J C K such that I fl (7(H) =: <Ji(H) ^ 0 
and dist(/, a(H)\l) = dg > 0. We say that <Ji(H) is an isolated spectral island, completely included 
in I. We notice that in this case there exist (j, N) e N* x N* such that (Ti(H) = U ■ 

We shall denote by J/ := {j + I,... ,j + N} the set of indices of the eigenvalues contained in the 
isolated spectral island (Ji{H) C I. 

Definition 1.6. Let us choose a smooth, closed and positively oriented simple contour C in the 
complex plain, which surrounds / and does not intersect the spectrum cr(H). We introduce the 
band projection and the band Hamiltonian: 

p :=Ej(H)= ^ J^(H-i)-^di, ni:=pn, := ^ j^l(H - (1.9) 

For j e N* we define the Bloch eigenprojections: p^ := ^ f.^ p(d)ddj where p(0) := 

\(j)j(9)){4'j(0)\ are 1-dimensional orthogonal projections in 

We shall learn from Corollary 12.21 that both operators p and have associated symbols Sp and 
resp. Sf, of class S~°°(X). In general, this is not the case for each individual pj (even for j G Jj). 
Also, the following statements hold true: 

[H,p,] = 0, HEi{H) = ^Dp(A,)p,. (1.10) 

j&Ji 


1.3 Our main results 

Assume that Hypothesis 11.21 is satisfied. Recall that a magnetic field is described by a closed 

N 

2-form B = ^ Bjk{x)dxj A dxk with Bjk defining an antisymmetric matrix that verifies the 
j,k=l 

equation dB = 0 (here d denotes the exterior derivative on differential forms). Let us point out 
that we make no periodicity or slow variation assumption on the magnetic field. 

Hypothesis 1.7. We shall only consider Bjk G BC°^(X) for any j < k. 

We shall consider the limit of weak magnetic helds, controlled by a small parameter e G [0, eo] 
(see Hypothesis 13.41 for more precise details). Roughly speaking we shall consider a family of 
magnetic fields {SejeGfo.eo] form B^^ := e where the components of 5° belong to BC°° (A) 

N 

uniformly with respect to e S [0,eo]. There exists an 1-form A := ^ Aj(x)dxj called a vector 

potential such that B = dA. This vector potential is not unique and any A' := A df with 
/ G C^(X) also verifies the equation B = dA'. It is well known that under our hypothesis on 
B one can always choose A such that Aj G C'“|(A). We shall fix a family of vector potentials 
{^e}ee(o,eo] having Components in a bounded subset of C^i (A) such that = dA°; then H,. = dA,. 
for Aj := eA°. 

The basic mathematical construction we use is the ^twisted pseudodifferential calculus’ (see 
(13.31) 1 introduced in [^[211I2S] (or [30] for its integral kernel version) and developped in [T711131133] . 
that associates a ’quantized operator’ Dp^(F) to any Hormander type symbol E G Sff(X). We 
shall use the shorthand notation Dp'^(F) := Dp^'(F). 

Recalling the results in let H'^ be the self-adjoint extension of Dp'^(h) with the domain 
given by the magnetic Sobolev space. Using the results in nms] or [2111117] we know that the 
resolvent set is stable for small variations of e. More precisely, given I as in Hypothesis 11.51 there 
exists eo > 0 small enough, such that cr(7L'^) fl/ 0 and dist(/, cr(iJ'')\/) > do/2 for any e G [0, eo]. 
Moreover, the Hausdorff distance between (j(H^) (~l I and cri{H) goes to zero with e. 

The main question we are concerned with is the following: if e G [0,eo] and cq > 0 is small 
enough, can we replace H by id" in (II.IHI) putting in the right hand side a power series in e 
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with leading term similarly defined but with instead of Op ? Our first result is contained in 
Theorem ll.QI and gives a partial answer to the above question. In order to state it, we shall need 
the following technical result: 

Proposition 1.8. For every 3 G p[H^) G denote the symbol of — 3 )“^ (as 

defined in m)- If l belongs to a compact subset K of C \ cr(H), there exists eo > 0 such that for 
e G [0, eo] we have that: 

1. K cC\(t{H^); 

2. the following development, convergent in the topology of the C*-norm || • \\w,e (induced from 
M{'H) as defined in (13.111) ) is true, uniformly with respect to G [0,eo] x K: 

= ^e'^rn{h-,e,:j), ro(/i; 6 , 3 ) = r„(h; 6 , 3 ) G 

raeN 


3. the map K B ^ i-p- r^{h) G '"(A’) is continuous for the Frechet topology on "^{X) 
uniformly in e G [0,eo]. 


Theorem 1.9. Under Hvvothesis M .‘A 11.51 and \3.4\ there exists eo > 0 small enough such that for 
any e G [0, eo] and for any n G N* we have that: 

1. H<^Ei{H'^) =Dp'^{Ssi) + e’^Op'^{vl)+e'^Dp%R}{h;n)) with: 

l<k<n—l 


(a) S^{x,0 ■■= 


It, E + y/2,e)(j)j{x - y/2,e) d9 


<jeJi 


dy, 


(b) vl := -(27rl) ^ Jchrk{h;e,i)di, 


(c) R}{h;n) :=-{2Tri) 3 ( T, e,i) ] d^. 


V k>r 


2 . there exists an orthogonal projection Pf ^ in L^(A’) such that 

(a) \\Ei{H^) - PfiJ < C„e", W[Ei{H^) - P/J|| < C„(/i)e", 

(b) ||[P%P/,J|| <C„(/r)e”, 

(c) = Op^(5p) + E e'=Op^(u|) + e^Qp^{Rj{p-n)) with 

l<k<n-l 


i. Sp{x,f) 


|S.| 




I, 


T* 


^ + y/2, d)(l)j{x - y/2, 0)^ 


de 


dy, 


a. ul ■.=-{2m) ^ J^rk{h;e,i)d^, 

Hi. R}{p; n) G 5'j“"™(A’)r uniformly for e G [0, eo]. 

Corollary 1.10. If the spectral island aj is the image of an N-fold degenerate isolated Bloch band 
A : T* —>■ R we have that 


H^Ei{H^)=Op^{\)Ei{H^)+ ^ e^Dp^{vl)+U^Dp^{R]{h-n)). 

l<k<n-l 


The second result is about the case in which the projection associated to the non-magnetic 
spectral island admits a basis of localized Wannier functions (see formula (I2.15|) for their defini¬ 
tion). Then we can associate to the band Hamiltonian a smooth N x N matrix-valued function 
{pjk{9)} defined on the dual torus T* (see Definition 12.71) : let us denote by fi its periodic extension 
to X* and by {(/ijfc)-),}.ygr its Fourier coefficients defining a sequence in ^^(r) having rapid decay. 
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We shall use the notation 


A^{x,y) := e (1.11) 

where the integral is taken along the oriented segment [x, y]. For every fixed pair 1 < j,k < N we 
consider the operator Op^Cfljk) defined in P(r) by the infinite matrix 

:= A (a, (1-12) 

The N X N matrix with entries Dpffor 1 < j,k < N defines an operator in P{T)^ that we 
denote by £)pp(/I). 

Theorem 1.11. Under Hvvothesis \1.2[ \1.5\ and \3.4\ if the isolated spectral band at zero mag¬ 
netic field admits an orthonormal basis consisting of composite Wannier functions (see Hypoth¬ 
esis \2.5\) . and if cq is small enough, then there exists an orthonormal magnetic localized basis 
{W^j}( 7 .j)GrxJ/ such that: 

1. = E E 

76riGJ/ 

2. sup sup(a: - 7 )™|>V^ (x)! < oo, V( 7 ,j) e F x Jj, 

7 erxGA' 

3. there exists a positive constant C < oo such that 

(a,/9)Grxr (j,k)GJiXJi 

Corollary 1.12. Under our assumptions, if cq is small enough, then the magnetic band Hamilto¬ 
nian H'^Ei{H^) is isospectral with Dpp(/i) up to an error of order e, i.e. the Hausdorff distance 
between their spectra is of order e. 

Our last result deals with the case of a slowly varying magnetic perturbation and makes the 
connection with the ‘usual’ magnetic Weyl calculus via minimal coupling. 

Definition 1.13. We say that a magnetic field Bfx) is slowly varying if it can be derived from 
a magnetic vector potential Afx) = A{ex) where A is a smooth vector potential having bounded 
derivatives of all strictly positive orders. In this case 

Bfx) = edA{ex), eG[0, cq]. 

Theorem 1.14. Besides the assumptions of Theorem \ 1.1 11 we consider a slowly varying magnetic 
field. Denote by Dp(/i'‘) the usual Weyl quantization of the matrix-valued symbol ■— 

Ujkif,—Afix)), acting onL‘^{X)^. Then the Hausdorff distance between the spectrum of H’^Ej{H^) 
and Dp{y,^) is of order e. 

2 The non-magnetic case 

2.1 Kernels and symbols of F-periodic operators 

Let us consider a general F-periodic bounded operator T G B(L^(A’)) (i.e. a bounded operator 
that commutes with all the translations {T.^}.ygr). Then Remark [T|3] implies that its Bloch-Floquet 
transform is a decomposable operator and we can write 

/*0 

i := '^pT'^p"^ = / ‘I{0)d0, T( 6 ») G V 6 » G T*. (2.1) 

Jt, 

We recall that the elements of are completely determined by their restriction to E C X 
and this restriction defines a unitary isomorphism .^g = L'^{E). Thus each fibre operator T(0) 


< Ce. 


WTG(X)) 
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has an associated distribution kernel T{9) € x X) that leaves invariant. Thus it verifies 

the relations: 

(T-a( 8 >r^)f( 6 ») = V 6 » e T, V(a,/3) e T^. (2.2) 

We shall always suppose that the map T* 9 0 i->- T(0) G ]B(^e) is bounded and measurable 
but much stronger hypothesis will be necessary in the sequel. We shall mainly be interested in 
situations in which the integral kernel of T(9) is of class BC°° {X x X^ and (12.21) implies that they 
satisfy the relation: 


T{9,x + a,y + f3) = e 




^>fi9,x,y), Vi9,x,y) gTxX^, V{a,l3) G T^, 


(2.3) 


Hypothesis 2.1. Suppose given a bounded T-periodic operator T such that in the Bloch-Floquet 
representation it defines a family of integral kernels T{9) of class BC°°{^X x X^ and such that the 
map T* 9 0 I—;> T{9) G BC°°[X x X) is bounded and measurable. 

Let us compute the integral kernel of T = Recall that every x £ X can be uniquely 

written as [x] + x with [x] G T and x £ E. For any / G S^{X) and x £ X we have 

{^f){x) = \E,\-^ [ d9{i{9)^rf){0,x) = \E,\-^ [ d9 [ dyf{9-,x,y){^rf){0,y) 
jt* Jt* Je 

= l^*r'E/ / dyT(0;i,y)e-*<^’M-^>/(y + 7) 

.^gpdT. Je 

= E/ + 7), 

d 0 f( 0 ;x,y)e-*<®’“>. 


7er ‘ 

t{a;x,y) := \E^\~^ 


(2.4) 




Thus T has the integral kernel: 

K%{x,y) := f([a;] - [y],x,y) := \E^\~^ 


f d0f(0;i,y)e-*<®'[“l-[^l>. (2.5) 

Je, 


Notice that this kernel may have some decay in the variable x — y depending on the regularity of 
the map T* 9 0 H> T{9) £ BC°° (^X x ff), which we shall assume to be (7°°. Indeed, by partial 
integration in ( 1 ^ we obtain that < x — y K%{x, y) is globally bounded for every N' > 1. 

Let us write down the symbol of the operator T 

S%{zX) = K<i{z + vj2, z — vl2)dv 

Jx 

= {2Tif J (^E^\-^ j d0f{0-z + {v/2),z-{v/2))^dv. ( 2 . 6 ) 

Proposition 2.2. //T is a bounded T-invariant operator satisfying Hvvothesis \2.1[ then its symbol 
3% belongs to S~°°{X). 

Proof. The result follows from ( 1 ^ and from the decay properties of K‘i{z + r/ 2 ,z — v/2) seen 
as a function of □ 

Proposition 2.3. Any function p £ C(T,) defines an operator Nip £ M{i^) given by multiplication 
with p{9) in each fiber space Then, denoting by p £ C'r(ft’*) the periodic extension of p to X* 
we have the equality: 

= Dp(p). 
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Proof. For any f G we have 




-1 


'T* 


p{0)de ^rf {x) = \E. 


/T* 


p(0)E' 


— i<6,'y> 


fix -i)\de 


7Gr 


= \E*GY. / /(x - 7) 


= ( 2 ^) 


— d 


X Jx* 


e*<«’^>p(^)/(x + z)(i6»(iz = (Dp(p)/)(x). 


□ 


2.2 The symbols of the ‘band operators’. 


Recall that the non-magnetic operator H has an isolated spectral island aj, see Hypothesis 11.51 
Let us list some properties of the band operators appearing in Definition 11.61 
Proposition 2.4. The following facts hold true: 

1. We have that [it, p] = [b,p] = 0. 

2. Using Provosition \2.2\ we may conclude that both operators p and Sj have associated symbols Sp 
and resp. Ss^ of class S~°°iW). 

3. Using the Bloch-Floquet transform (see subsection \1.2\) we can write the band operators as direct 
integrals of some fibre operators given in terms of the Bloch eigenvalues and eigenvectors: 



and we have the following formulae for their integral kernels and symbols: 

Kf,ix,y)= E {\E*G [ Xkid)<(kix,0)(j)k{y,0)d9 


kGJi 


/T* 


i^p(a;,y) = E (t>kix,0)(j>kiy,0)de^ 

Ss^izX)= E (i‘^^ f\E*G f f Xki0)Mz + v/2,9)Mz - v/2,9) do) dv, 

k&Ji X dx dT. / 


SpizX)= E f f Mz + v/ 2 ,e)Mz-y/ 2 , 9 )d 9 )dv. 


(2.7) 

( 2 . 8 ) 

(2.9) 

( 2 . 10 ) 

( 2 . 11 ) 

( 2 . 12 ) 


We can define for any j S J/ the Bloch orthogonal projections pj(0) := \(l)ji9)){(j)j{9)\ and 
pj := ^ pj{9)d9^ '^r- Let us notice that in general these one dimensional projections have 

much less regularity in the variable 0 € T* than the finite sums above, and thus these objects are 
not very useful without some stronger hypothesis. 

In the case of an isolated spectral band we have that p = X) Pj ^bo 

jeJi 

HEi{H) = E^P(Aj)p,i?/(ii), (2.13) 

j&Ji 

e-^*^Ei{H) = E^pE^'^OPj^H^)- (2-14) 

j&Ji 
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2.3 The composite Wannier basis. 

Hypothesis 2.5. Under our Hypothesis 11.51 of an isolated spectral band, we suppose also that 
the subspace Tg := p{9)^e (see Definition 11.61) has an orthonormal basis tjjj{9) G H BC°°{X) 
indexed by j G J/ and such that the map T* 9 d >->• G BC°°{X) is smooth for any j G Ji- 

We do not suppose that the vectors of the basis are eigenvectors of H{9); but we notice 

that they are finite linear combinations of such eigenvectors and thus smooth functions ot x G X. 

We have that 'ipj{9) G p{9)^e C V(^H{9)Y This smooth basis allows us to construct a so-called 
composite Wannier basis for the range of p = Ei{H). 

Proposition 2.6. Define the functions 

Wj := e y{X), Vj G Jj; Wj{x) = \E^\-^ [ (2.15) 

that are orthonormal in Lfi{X). Then the translated functions W-y.j := T-^Wj with (7, j) G T x J/ 
form an orthonormal basis for pTL. Moreover, < x >*” vjj{x) lies in L°°{X) for all m > 0. 

Proof. We note that ('^rW 7 j) (i, 0) = {x,9) and the proposition follows from the fact 

that {'4!j{9y}j^jj is an orthonormal system in for each 9 G T*. 

We assumed that each ij}j{x,9) is smooth in 9 which implies that < x >™ Wj{x) lies in Lf{X) 
for all m > 0. Moreover, because the integral kernel of p obeys the estimate \Kp{x,y)\ < Cm < 
X — y valid for every m > 0 , it follows that < x >™ Wj{x) lies in L°°{X) for all m > 0. □ 


Definition 2.7. 

1. We define the effective Hamiltonian associated to cr/ to be the N x N matrix valued map 

T* 9 d fijki9) ■■= {il}j{9),H{9)il}k{9))^g G C, (j, fc) G Ji X Jj. 

2. For any {j, k) G Ji x Jj we define the rank-one operators: TTjk{9) ■= \ipk{9)){ipj{9)\ and their 

associated F-invariant operators in L'^{X) given by IVjk '■= We shall denote by 

Pjk the integral kernel of the operator H^-fc. 

Remark 2.8. We have the following properties: 

1- {'^jk}(j^k)eJf i® ^ family of bounded operators on that satisfy the relations: 

^jk ~ ^kj ■ ^ J I '^jk'^pq — djqTTpk- G Jj. 

2. Similar properties are also valid for the family {Iljk}(^j^k)eJ] ■ Moreover, for every pair 
(j, fc) G Jf, there exists a symbol pjk G S~^{X) such that Hj^ = Dp{pjk) and we have 
the explicit formula 

Pjk{z,C) = {2TTf\E.,\~'^ f [ ipk{z + iv/2), 9)iIjj{z - {v/2),9)d9dv, V(z,C)gS. 

3 . p = andpHp = /® X) yjk{9)TTkj{9) d9. 

jeJi * (3,k)ejf 

We also list without proof a few properties of the band Hamiltonian seen as a matrix in the 
Wannier basis: 

Proposition 2.9. The following equalities are true (here Jljk is the periodic extension of fXjk): 

p = N~^ Op{pjk)Op(j)kj)-, b = pi7p = N~^ Y ^PiT3k)i^p{pkj); (2.16) 

U,k)eJf (i.fc)eJf 
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(y^OL,j 5 n/j7j,W^ 

(W„,, , HWp,k)L^^x) = l^*r' I =: 

(W,,, , = \EA~" 

3 Adding a weak magnetic field 

3.1 Brief recall of the magnetic Weyl calculus. 

We shall very briefly recall some of the main definitions and results concerning a completely gauge 
covariant version of the ‘minimal coupling’ procedure developed in p51 [T71 [T51150] . 

Given a magnetic field B and a choice of a vector potential A for it and considering the 
fundamental set of dynamical observables given by the 'minimal coupling' hypothesis: 

{Qi,...,Qd}, {Qjf){x) := Xjf{x), yf € y(X) (3.1) 

{nf,...,n^}, {uff){x):={-id,f){x)-A,{x)f{x),yf€A^{x), (3.2) 

in m we considered the twisted Weyl system defined by the unitary groups associated to the 
self-adjoint extensions of the above operators. This procedure allows us to define a ‘twisted pseu¬ 
dodifferential calculus' (introduced in [531 [Ml US] and developped in [T71[TH1[M]) that associates 
to any Hormander type symbol F G SJf(X) the following operator in Lf{X) (for all u G dX{X) 
and X G X)\ 

{£)p'^{F)u){x) := {2Tr)~^ J j e~"^F{^^-^,^)u{y) dfdy. (3.3) 

Two important results in [23j state that two vector potentials that are gauge equivalent define 
two unitarily equivalent functional calculi and that the application defined above extends to 
a linear and topological isomorphism between JX'{E) and (X); JX' {X)) . At this point we can 
make the connection with the ‘twisted integral kernels' formalism in [30] . where for any integral 
kernel K G dX'{X x X) one associates a twisted integral kernel (see (11.91) ') 

K^{x,y) := A^(x,y)K{x,y). (3.4) 

For any integral kernel K G ^'(X x X) let us denote by IntAT its corresponding linear operator 
on A^(X) (i.e. (FntFu)(x) = K{x,y)u{y) dy). Let us recall the usual Weyl calculus, that 
we shall denote by Dp and the linear bijection 2IJ : y'{^) —> dX'{X x X) associated to it by 
Dp(F) =Xnt(2rrF): 


(2BF)(x,y) := {27r)-^ j^^F<^’-y>F{^,^)d^. (3.5) 

Then we have the equality 

Dp^(F) = Xnt(A^2IIF). (3.6) 

Remark 3.1. We notice that under our Hypothesis 11.71 and with the notation ()1.9I) : 


A^{x,z)A^{z,y)A^{y,x) = exp / E V =: n^{x,y,z); 

I J<x,y,z>) 

the above integral is taken on the positively oriented triangle < x,y, z >. We have the estimation: 

\^^{x,y,z) - l| < CIlHlIooKy -a;) A (z - x)\. 
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Using Theorem 4.1 in [m, under Hvpothesis ll.71 for a symbol h G S'^{X)y- verifying Hypothesis 
11.21 the operator Op‘^{h) for any A with components of class C^^{X) has a closure in L‘^{X) 
that is self-adjoint on a domain "H™ (a ‘magnetic Sobolev space’) and it is lower semibounded. 
Thus we can define its resolvent ~ i)~^ for any 3 ^ a{H^) and Theorem 6.5 in [T5] states 
that it exists a well-defined symbol r^(h) G S^'^^X) such that 

= Dp^(rf(h)). 

Let us also recall from [53] that the operator composition of the operators D\)^{F) and Dp'^(G') 
induces a twisted Moyal product depending only on the magnetic field B: 

{FfG){X) := dU ^F{X - Y)G{X - Z) (3.7) 

where we use the notation of the type X = {x,^),Y = {y,r]) etc. and we have denoted by T{x, y, z) 
the triangle in X of vertices x — y — z„ x + y — z, x — y + z and by T{x, y, z) the triangle in X 
of vertices x — y + z, y — z + x,z — x + y. For any symbol F we denote by Fg its inverse with 
respect to the magnetic Moyal product, if it exists. 

We shall often use the symbol < ^ >™ for m > 0 and we shall also need to consider its magnetic 
Moyal inverse. For that we use the arguments in Section 2.1 of [5^ and conclude that for a > 0 
large enough the symbol 5m{x,^) ■=< ^ >"* +a, with m > 0 has an inverse for the magnetic 
Moyal product and we shall use the shorthand notation 5 ^^ instead of {sm) g for this symbol 
with some fixed large enough am > 0 . 

For the completeness of our arguments we give the proof of a simplified version of Proposition 
8.1 in m that will be important in our arguments. 

Proposition 3.2. Suppose we are given (j) G Sp(S), G S'p(S) and 0 G BC°°{X;G^^{X'^)) 
(the bounded smooth functions on X with values in the space of smooth functions on X^ with 
polynomial growth together with their derivatives). Then 

£(6»; (j), fj){X) := J dY dZe-'^^’^°^^'^'>e{x, y, z)(l){X - Y)'4 j{X - Z) 

defines a symbol of class and the mapping 

5™(S) X s^iE) 9 (</),v-) ^ e 5p"^+^(s) 

is continuous. 


Proof. By a simple change of variables we can write 

£( 6 »; (j), i()(X) = J dY x-y,x- z)(j){Y)i({Z). 

For any natural numbers fVi, iV 2 , Mi, M 2 we have the identity: 


1 + 2|^-CP 


-2ic!°(x-Y,x-z) _ ( 1 ~ ~ C)i dy) \ ^ / 1 + i{{£, — v)fdz) 


l + 2|?-r7|2 


l + i{{x- z),d^) \^^ /^ l-^((a;-j/),^c) ^”' - 2 ia°ix-Y,x-z) 


Ni 


Ml 


(3.8) 


l + 2\x — z\'^ J \ l + 2\x — y\‘^ 

If (j) and Ip are test functions, after integration by parts we obtain the estimate: 

|£( 6 l; (l>,i^)iX)\ <c(^jjY<£,-y >-^^< 77 >™< x-y 


dZ <^-C >-^=< C>P<x-z> 


r2{Ni,N2) — M2 


<G' <f> 


m+p 
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where we choose > |m| +n, N 2 > \p\+n, Mi > ri{Ni,N 2 ) + n and M 2 > r 2 {Ni,N 2 )+n, with 
rj(iVi, 7 V 2 ) the powers dominating d^^dy‘^0{x, y — x, z — x). Now let us compute the ^-derivative 
of 


(X) =-y dY J dZ + d(;.)e 9{x,x - y,x - z)(l){Y)^{Z) 

= X dY jjZe-^^^°^^-^'^-^^0{x,x- y,x-z) [{dy^cj,) (F)] V'(^) 

+ J dY J^dZe-^^'^°^^-^’^-^^B{x,x - y,x - z)ct){Y) [(S^.V') {Z)] 

= £(0;(%(^),^)(X) + £(0;<^, (%V^))W. 

Considering the x-derivative we obtain in a similar way that 

where 0{x,x — y,x — z) := dxjd{x,x — y,x — z). These two formulas allow us to control all the 
seminorms in the corresponding Hdrmander symbol spaces. □ 


Fix x,y,z G and define the following objects: 


Dfk{x,y,z) := f ds f dtBjk{x + {l-2s)y+{2t-l)z), 
Jo Jo 


(3.9) 


FB{x,y,z) :=^ f B = {D^ {x,y,z) z, y) := '^yjZkDf^{x,y,z). (3.10) 

4 JT{x,y,z) 

Corollary 3.3. For every (m, m') € K x R and p G [0,1], and for any magnetic field B satisfying 
Hypothesis \1.1\ we have that the map x (M'^) 9 (a, &) ,a^^b G 5'™+™ (R'^) is 

bilinear and continuous. 


Proof. We use the formula of the magnetic composition (EID and apply ProDOsition l3.2l bv replac¬ 
ing e{x, y, z) with . □ 


Using Theorem 3.1 and Remark 3.2 in m we deduce that for any $ G =^(S) the operator 
Dp'^(<i)) defines a bounded operator on L^(A’). This allows us to define on =5^(S) a C*-norm that 
only depends on B: 

mW.B ■■= ||Op-^(‘i>)||B(L=(A')), Vd> G S^{E). (3.11) 

Using the above cited results in m we can extend the above C*-norm to the space Sq(Y). 
Moreover, it is also proved in m that this norm is bounded by a specific norm from the family 
defining the Frechet topology of Sq(Y), that depends only on the dimension d = dimRfh. 

3.2 Weak magnetic fields 

In our paper we are interested in weak magnetic fields, that we shall control by a small parameter 
e G [0, eo]- 

Hypothesis 3.4. We shall consider a family of magnetic fields {Rejjgjo.eo] form B^ := e 

with Bg having components in BC°° (ft) uniformly with respect to e G [0, cq]. In order to simplify 
the notation, when dealing with weak magnetic fields, the indexes (or the exponents) or 
shall be replaced by e. Similarly, we shall use the notation || • \\w,e instead of || • 
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Proposition 3.5. For e G [0,eo] there exists a continuous application r^ : x (X) — >• 

such that: 

af6 = afb + er,{a,b), V{a,b) G 5'™(K'') x sf(R^), (3.12) 

Proof. In (I3.7|) we use the following identity: 

g-4iF, ^ ^ f 

Jo 

in order to obtain (I3.12|) with 

[r<:(a,6)](X) = i(af& - afb) (3.14) 

= {D^{x,y,z)z,y)a{X -Y)biX - Z)dYdZ. 


All the components of the matrix belong to BC°°[X) uniformly with respect to e G [0,eo] for 
any j < k. Integrating by parts we obtain: 


[r,(a,6)](X) = J 

d 

^ H -Y) {d^^^b){X - Z). 

j,k=l 

The proof can be completed by using Proposition 13.21 where 


9{x,y,z) 


pAitF^(x,y,z)^^ 


D%{x,y,z). 


(3.15) 


□ 


Remark 3.6. If we replace (13.131) with the A^’th order Taylor expansion of the exponential, we 
obtain that for any iV G N*: 

afb = afb + ^ e^cf\a,h) + e^ pf\a,b), (3.16) 

l<fc<Af-l 

with cf\a,b) G 5'™+'" and pf^\a,b) G 5'™+'" ~'^^^[X) uniformly in e G [0,eo]- 

Associated to the series development of the symbol r'^ih) given in Proposition 11.81 we shall 
also use the notations 

rln(.h)-.= ^ eVfc(h;e,3) G5r™(A)r; ^ eVfc(/i; e, 3 ) G Wr- (3.17) 

0</c<n n+l<fc 

Remark 3.7. The remainder G 5')“™(A)r has the following properties: 

1- where G S')“’”(A)r uniformly in e G [0,eo) for some small enough 

Co > 0; 

2. hfr^j^^n = e'^^^hfr^fi ^ ( — re(h, r°)^'*) j and noticing that+ 

yfc—n+1 j 

erfh,r^) we conclude that < (76"+^ for some C > 0 independent of e G [0,eo) 

for some small enough cq > 0. 
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3.2.1 Proof of Proposition ITTSl 

In this paragraph we shall use the results above in order to prove Proposition [T^ The first point 
clearly follows from the spectral stability results proved in mi, as briefly recalled in subsection 
11.31 before the statement of Proposition 11.81 

For the last two statements we start from the continuity of the application 


C\<j(R^) 3 3 ^r^(h)eSr"^(A:) 

for the II ■ II w,e-topology, a consequence of basic spectral theory for self-adjoint operators. In order 
to obtain a control in the Frechet topology on Si"‘(A^) we recall some results from [T5]. Let us 
recall the symbol Sm introduced in Subsection 13.11 lust before Proposition 13.21 and the space of 
“linear” symbols: 

vxes, lx(Y) :=a°(X,Y). 

For any X G 5 we can define the operators 

a^xbP] := Ixi'-'tp - e y'{E). 

Then Theorem 5.2 in [T5] states that the Frechet topology on any space (for any m G R.) 

may be also defined by the following family of seminorms: 




s„+gf • • • ac);(JV’]) 


G R+ 

W,e 


(3.18) 


indexed by a pair of natural numbers (p,q) G N x N and by two families of points {ui, ..., Up} C X 
and {pi, ..., pq} C X*. A simple computation shows that for any e G [0, eg] and any 3 ^ o'(iL'^) 


= -rl{hWad^x{hWrl{h). (3.19) 

Using the resolvent equation: 

rlih) = rl{h) + {i - 3 )rl{h)rl{h), (3.20) 

and Propositions 3.6 and 3.7 from [18] we easily prove that the applications: 

K e S^{X) (3.21) 

are well defined, bounded and uniformly continuous for the norm || • \\w,t any e G [0, eg]. 

The second point follows by noticing that the result in Lemma 13.51 implies the equality 

1 = ih- 3 )fr°{h) = (/i- 3 )fr°(h)-Fer,(/i,r°(/i)) (3.22) 

with the family {r'e(/i,r°(/i))}gg[o_Eg] being a bounded subset in S~^{X). We conclude that for 
some eg > 0 small enough, l-l-ere(h, r°(/i)) defines an invertible magnetic operator for any e G [0, eg] 
and its inverse will have a symbol s'^( 3 ) given as the limit of the following norm convergent series: 

5 ^( 3 ) := E( - G S°{X). (3.23) 

rtGN 


This clearly gives us the development in point (2) of the Theorem with 

r„(/i;e,3) := (-l)V°(h)f (r,(h,r-°(h)))*'” G (3.24) 

In order to control the uniformity with respect to e G [0, eg] of the continuity of the application 
in (|3.21|) let us notice that 

rlih) - rpih) = { 3 ' - 3)rl{h)i^rl,{h) (3.25) 

and that for any 3 G AT the family of symbols {?’|(^)}eG[o.eo] > eo > 0 small enough, is a bounded 
set in S~"^{X) due to point (2). 
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3.3 Proof of Theorem 11.91 

We shall consider an isolated spectral band as in Hypothesis II.51 with a hand projection and a band 
Hamiltonian as in Definition 11.61 they describe the exact dynamics in the given ‘energy window’ 
in the absence of a magnetic field. Now suppose that a magnetic field satisfying Hypothesis 13.41 
(depending on the parameter e € [0, eo]) is switched on. 

We introduce the following simplified notations: 


p" := ■= H^Ei{H^). (3.26) 

We have that , p*’] =0 and 

P' = h" = ^ (3-27) 

with C the smooth closed curve in C defined just before Remark 1 1.91 for H replaced with 77“ with 
e e [0,eo] and cq > 0 small enough. 

Proposition 3.8. Under Hvvothesis \1.5\ and \3.4\ and using (13.171) . we have that 

ij“=Op^(^^)+ ^ eW(^fe)+eW(7?!(/i;n)), (3.28) 

l<fc<n—1 


with 


Ss^izX) = 




>^kid)4’kiz + v/2, 0)(j)k{z - v/2, 0) d9 dv 


(3.29) 


and 


■=-{2m) i irk{h]e,i)di, R}{h;n) 



Proof. We use (13.271) and the development of r^{h) given in the second point of Proposition 11.81 
Formula (I3.28|l follows from (11.911 and (12.111) . □ 


This Proposition proves point 1 of our Theorem 11.91 For the proof of the second point let us 
consider the projection p“ =: Dp“(S'pe) and try to approximate it by the magnetic quantization of 
the ‘free symbol’ S'p. A technical difficulty comes from the fact that now Dp“(S'p) is no longer a 
projection (it is idempotent only modulo an error of order e!). 

Proposition 3.9. Under Hvvothesis \ 1. 5\ and \3.4\ with the above notations, for any e € [0, cq] with 
eo > 0 small enough and for any iV S N* there exists an orthogonal projection p^ such that: 

Pn = J^rm{h;e,d)di) + e^Xf,j{h), with \\Xf,j{h)\\ <Cn < oo for any 

0<m<N—l 

e S [0, Co]. 

2. ||p“ — p^ll < CN{h)e^ for some CN{h) < oo. 

3. ||.5'’ — 77“p)^^|| < CN{h)e^, ||.Q“ — < CN{h)e^ for some Cnih) < oo. 

I II , P^] II < Cnc^, II [H^ , P^] II < Ve e [0, eo) 

5. 3RN{h) e S-°°{X) such that H'^p% = + e^Dp"(i?Ar(/i)). 

Proof. Using once again (13.2711 and the development of r|(/i) given in the second point of Propo¬ 
sition EM we can write 


5 


p' 




s, 


pym ■= ^ / 'rm{h,e,i)di. 


(3.30) 
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We define for any N G'N*: 


qiN) _ 

P' ■“ 27rz 


-J^rl^ih)di, zW:=V-5, 


(JV) 

p' 


^ := ( 1 / 2 ) ^ Zi^'> := ( 1 / 2 ) = Sp. - k 


(AT) 
P' ■ 


(3.31) 

(3.32) 


Then 


i G S^^{X), Dp^(zW) < C^e^+I. 


We define 


:= Dp^(5W), 

and notice that it is a self-adjoint operator and 


(3.33) 


(°w^r2_°(^) ^ = zW-(zWr5pe+5p.fzW) + (zW) 


f2 


implying the estimate ||(q^)^—q^|| < Ce^+^. Following the procedure in (Proposition 
3) we notice that this estimation implies that there exists ei G (0,eo] small enough such that 
tT(q'^) = IqU h where Iq C [—6-^+^, and /i C [1 — 1 -I- for any e G [0, ei]. Thus 

1 — ^ cr(q^) for any e G [0, ei]. We can thus find a smooth contour Ci C C such that 

[1 — 1 -|- be in in the interior region and [—6^“'"^, in its exterior region, and define 


(q^-i - 3 ) 


(3.34) 


By definition, p^ is a self-adjoint projection that commutes with q^_i; moreover it is equal to 
^^/i(qAr_i)- This means that 

Pn ~ Pat-i = Pat ~ Pw-i(pAr + (1 — Pw)) = (1 ~ PAr-i)pAr + qw_i(]l — P^) 

satisfies the estimation ||p^ — q^_i|| < Ce^. 

From (13.341) it follows that 


= -r. 

2711 


= "ii, Fp- -i):r(sr"'’ -v)f(Sp- 

= *'• + 2^ i, Fi"'" -3)J‘ ( 2 ^/pV,»-.W<ij') f (Sp. 

= + Si i F'- ->):«■ (i /yv.^-.wpis') f («<?-« -3):<i3. 

Using the second point in Remark 13.71 we easily obtain the following stronger result: 


3eo > OVIV G N* 3C]s[ < 00 


hf(%-5pe) < VeG[0,eo). (3.35) 


W,e 


We also note that [h , Sp, = htl‘^(5pe — 5'pe) — (Spe — Spe^'i'^h. From the bounds: 


5p^ - 5p 


(Af-l) 


W,e 


= ||p^ - q^-i|| < Ce 


N 


(3.36) 


we conclude that the following estimates hold: 


hr{Sp. - Sp.) < Cnc^, Ve G [0,eo). 

^ W,e 


(3.37) 
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Since 


we can write: 


||q^ - P'll < ^ jjr<^i,N\\w,e\dd\ < Ce 


N+l 




< Ce^, 


W,e 




W,e 


< Ce^. 


Using the magnetic symbolic calculus m we end the proof of the Proposition. 


(3.38) 

(3.39) 

□ 


Corollary 3.10. Under Hvvothesis I f . .51 and \3.4[ with the above notations, for any TV £ N* there 
exists eo > 0 such that for any e G [0, eo] we have: 


e-**" Ei{H^) = + e^D\'%{h) 

—-e 

with ||91'^(h)|| < CArltj for some Cn < oo and any e £ [0,eo]. 


(3.40) 


Proof. In order to prove (13.401) let us define the application: 

RBt^ 4'(t) := £ M{L^{X)) 

which due to (13.391) is norm differentiable for alH £ K and has the derivative: 

vi/(t) = -ie-^*^‘ e**^' £ MiL^X)). 

From (13.391) we get \\H^p% - Pat-^*'118(^2< Cnc^ and also 

<CN{h)\t\ 



< 

f 

i'(s) 

ds 


B(LUX)) 

Jo 


mux)) 


,N 


If we define 

^t) := 

and an argument similar to the one above, using now the estimation for , p^] allows us to 
conclude that 

<C^(/i)|t|e'^. 




mux)) 


Hence: 




MiL^X)) 


<CNih)\t\e^, 


which together with ||p^ — P'^H < C]s[{h)e" it ends the proof of the Proposition. 

3.4 Proof of Theorem 11.111 


□ 


Having the composite Wannier basis at our disposal (see (12.151) 1. we shall construct a ‘magnetic’ 
Wannier basis for the magnetic band projection as in [5], based on the strategy in [27) . 

The integral kernel corresponding to Ei{H) may be written as 

K{p){x,y) = ^ -j)wj{y -j) 

76 r \jeJi 

where the series is absolutely convergent due to the spatial localization of the Wannier functions 
(see Proposition l2.6l) . 

Definition 3.11. We introduce the following modified band projection kernel 


Ke{p){x,y) := '^\'^ujt^{x,y)wj{x--f)wj{y--f)\, ujt^{x,y) := (3.41) 

7er \jeJi J 

and the following bounded self-adjoint operators on L‘^{X): 
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1. := Int[A^if(p)] = Dp^ (5p), 

2. q" :=Int[A"A:,(p)]. 

We remark that 


K^{x,y)K^{p){x,y) = ^ '^A^{x,-f)wj{x-j)A^{y,'y)wj{y-'y) . (3.42) 

7Gr \j€Ji J 

Definition 3.12. Let us define the following modified Wannier functions (similar to those used 
in [I1I27115] for the constant magnetic field case) 

yV^^jix) := A%x,j)wj{x -j). 

Lemma 3.13. For any m S N there exists a eonstant Cm < oo sueh that 

sup <x-y>'^ \K^{p){x,y) - K{p){x,y)\ < Cmi- 

{x,y)eXxX 

Proof. First we fix some 7 S F and notice that for any {Ni,N 2 ) € there exists a constant 
C < 00 such that 

sup < a; - 7 y - 7 |a;^(a;,?/) - l| |uij(a; - 7 )||uij(?/- 7)1 < Ce. (3.43) 

(x,y)eXxX 

This allows us to conclude that for any m G N and any p G N 

< X-y >'^ lK,(p)(x,y) - K(p)(x,y)l < eCm,p sup V < a; - 7 >"^’< y - 7 

{x,y)eXxX 

which ends the proof. □ 


Proposition 3.14. There exists some C > 0 such that 

llq'^ ~ P*^!! < eC, Ve G [0, Co]. 

Proof. If we replace A^ = 1 in Proposition 13.91 (1) and (2), and using Definition 13 .11 1 we obtain 

llp^-q^ll <Ce. (3.44) 

From our Lemma 13.131 we get that — q*^!! < eC and the triangle inequality finishes the 
proof. □ 


Since A'^{Q,j)T-.y is unitary, for every fixed 7 G F the set {Wf is orthonormal, but this 

is no longer true when considering pairs {yVa,jj^p,k) a ^ /?. We shall apply an orthogonal- 
ization procedure in order to obtain an orthonormal basis for the range of q*^. 


Lemma 3.15. With the above notations, by writing (yVa k/ 2 ( ) ~ ^ap^jk + ^\^'ap\jk’ 
then for any m G N there exists Cm > 0 such that for any e G [0, eo] 


/^r 


<c„ 
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Proof. 


= j J^''{a,x)wj{x - a) h.^{x,l3)wk{x - I5)dx = 5ap5jk+ (3.45) 


lx 


+ A'^(a,/3) / PP{x,a, I3)wj{x — a)wk{x — (3) dx 


where for any m G N 


-/sr 


/ fl'^{x,a, I3)wj{x — a)wk{x — (3) dx 

Jx 


< 


<eCe/ < X — a > 


m+1 


/X 


'j{x — a)\ < X — j3 \wk{x — I3)\ dx < e 


for some C!^ ,, > 0 uniformly bounded for e £ [ 0 , eg]. 


(3.46) 

□ 


We shall follow replacing the one-dimensional situation considered there with our N- 
dimensional situation. Our family of modified Wannier functions can be seen as elements in 
Z“(r;L2(;t’;C^)). 

Lemma 15.151 can be restated as = 1-1- eX*^, where ^ := , ^ and 

p £ B(C'^) for {a, 13) £ L x L. Both families and X*^ define bounded operators both on 

Z^(r;C'^) and on uniformly for e £ [0,eo]. Thus, as bounded operators on Z^(r;C'^), 

1/2 

for eg > 0 small enough, we can define (G*^) as norm-convergent power series in e £ [0, eg]. 
Lemma 3.16. We can write 

(G‘=)“^/^ = 1 -h eY’= (3.47) 

with ^ £ B(C^) such that for any to £ N there exists Cm > 0 such that for any e £ [0, eg]; 

1/2 

Thus (G*^) can be extended to a bounded operator on Z°°(r;C'^), uniformly for e £ [0,eg]. 

Proof. We just have to use the norm convergent series expansion of the operator (G*^) and 
notice that for each term of order fc > 2 we can write 

<0-/3 >”"< < a - 7i >”"< 7i - 72 >""... < 7fc-2 - 7fe-i >”"< 7fc-i - /3 >"" • 


□ 


We define the family 

W := (G'^)"^^^W'^ £ L°“(r; [l2(A’)]^). (3.48) 


O O _- 

Its components {W^j}('y,j)erxJi form an orthonormal family and |]VV^j — < Ce uniformly 

in ( 7 ,j) £ r X J/ and e £ [0,eg]. Moreover, for any 7 £ T and to £ N there exists Cm > 0 such 


that 


sup < a; — 7 >™ 
xex 


W^j{x) 


WfAx) 


<C„ 


(3.49) 


uniformly in ( 7 , j) £ T x J/ and e £ [0,eg]. 

Definition 3.17. Let us define the orthogonal projection 


Proposition 3.18. 
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There exists C > 0 such that ||p*^ — p'^|| < eC for any e € [0, eo]. 


• For any e € [0,eo] the operator ■= t — (p*^ — p*^) [p'^p'^ + (l — P*^) (l — P*^)] ts 

a unitary operator satisfying the intertwining property '^'^p'^ = p'^'^'^. Thus the family 

O 

j)&rxJi t-s an orthonormal basis of which is an invariant subspace of FT^. 

• We have the following norm-convergent power series expansion: 

= 1+ Y. e™T„(e), 

m^N* 


where all the operators Tm(e) are uniformly bounded in e G [0,eo]- 

Proof. We use (13.4911 in order to estimate the operator norm of p*^ — and show that it goes to 
zero like e. Then we use Proposition 13.141 and the Sz. Nagy formula of the unitary intertwining 
operator (see Remark II.4.4 in [10]). □ 


We are now ready to define the ‘magnetic Wannier’ functions as the following family indexed 
by ( 7 , j) G r X J/: 

3.4.1 Concluding the proof of Theorem 11.111 

The first two points of the Theorem are direct consequences of the definitions and arguments 
above from which we also conclude that for any m G N there exists a finite positive constant Cm 
such that 


sup < a — fi >'' 
(a,/3)erxr 




L^iX) 


C: Cm£ 


which is independent of the indices (a, /3) G T x T and (j, k) G Jj x Jj. Now we compute: 




L^(X) 


Wj{x — a)A'^{x,a)A'^{x,x')A'‘{fi,x')Ksj{x,x')wk{x' — fi)dxdx' 


a: Jx 


= A^(a,/3) 


n'^(a, X, x')VF{l3, X, x')Kfj(x, x')wj{x — a)wk{x' — j5)dx dx' 


X Jx 


— A {CX, + {Xj j.)al3 — iDp-p{p.jk)a,l3 F {Xj i.)ai3 


where in the last equality we used (11.121) . Using Remark 13.11 and the fast decay of the Wannier 
functions {wj{x)}j^jj, we obtain for any to G N the estimate 

sup <a-/l>’" |(X;fc)„; 3 | < (3.50) 

(a,/3)erxr 


uniformly in (j, fc) G Ji x Jj and e G [0, eo]. This implies that for every to > 1 there exists a 
constant Cm > 0 such that 




< CmS < a — B >~ 


(3.51) 


In particular, this estimate proves the third point of Theorem 11.111 and the spectral result given 
in Corollary II. 121 
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3.4.2 An alternative form for the effective Hamiltonian. 


We shall cast now the result in Theorem 11.111 in a form that generalizes formulas (I2.16P to the 
case of a small, smooth magnetic field. 

We consider the Hilbert space 'Hn '■= C^, with N gN* from Hypothesis 11.51 

We denote by Iat the identity operator on C^. 

We notice that due to the decay properties of the ’magnetic’ Wannier functions introduced 
above, using the Cotlar-Stein Lemma allows us to prove that for any pair (j, fc) G J/ x J/ the 
following series are strongly convergent in B(7^) defining bounded operators (with good estimations 
on the norm): 

( 3 - 52 ) 

It is easy to see that H^-^ are orthogonal projections for any j G J/ and we have the relations: 

[n^\]*=n^,., = (3.53) 

Let us define the operators H*^ G M{'Hn) given by the B(H)-valued N x N matrix with entries H^^ 
and Op’^Cfl) given by the B('H)-valued N x N matrix with entries for {j, k) G Jj x Jj. 

Proposition 3.19. Under the Hypothesis of Theorem M.lll there exists some constant C G K+ 
such that for some small enough eg > 0, we have for any e G [0,eo] that: 


Sd" ® Iat 




< C7e. 

B(Wn) 


This Proposition follows imediately from our Theorem 11.111 and the following statement. 
Proposition 3.20. For an isolated spectral band / C M that also satisfies Hvvothesis 1 2. .51 and for 
the linearly independent system {Wq j}(a,j)erx J/ defined above, the following relation holds: 


for any {a,j3) G T x T and any {l,m) G Jj x Jj, with some positive constant that is uniformly 
bounded for e G [0, eg]. 


Proof. It is evidently enough to compute the following scalar products: 





/ A^O, a :)7 {x - a) A^a:, x + 7 ) P)T-i3f>rn) (x) dx 

7Gr 

J ?c 


/ A^{a,x)f)i{x - a)A'^(x,x + 7 )A 7 x + 'y,P)tprnix - P + ^)dx 

7Gr 

J X 


+ X !(^)-7 / [f2"(a,a:,x + 7 )H'’(a;,a; + 7,/3) - l]tjji{x - - f3 + j) dx. 

7er 

In order to estimate the last contribution above we notice that 

X!|(^)-7 / \^''{o:,x,x+ 'y)n''{x,x+ 'y,P) - l\\tjji{x - a)\\tfjm{x - p+ j)\dx < 

<eCe'^\{]Mfn) < 7 >^ / < X - a > l-ipi^x - a)\ < x - P + j > l-ipmix - P + j)\dx < eC',,. 

7er' 

□ 
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3.5 Proof of Theorem 11.141 

Let us recall that when composite Wannier functions exist, the non-magnetic band Hamiltonian 
HEi{H) is unitarily equivalent with a matrix in P{T)^ (see Proposition 12.91) given by 
where a, /3 GT and 1 < j, k < N. 

We now investigate what happens when the magnetic field perturbation is not just globally 
small but it also has slow variation, i.e. it is generated by a vector potential Ae(x) = A(ex), where 
A has bounded hrst order derivatives. From Corollary 11.121 we know that the magnetic matrix 
is (up to an error of order e in the norm topology) unitarily equivalent 

/ L^{X) 

with: 


exp 




(3.54) 


We compute: 


f A^= 'Y] {l3-a)j [ dt Aj[e{a +l3)/2 + te{l3 - a)) 

l<j<d J-1/2 

^ A AA 

= {A^{{a +/3)I2),I3 - a) + e^{l3 - a)j{/3 - a)k ds tdt{dkAj) {est{l3 - a)). 

Jo J-l /2 


rl/2 

3,k •''J 0-1/2 

Taking into account the hypothesis on the derivatives of the vector potential we obtain: 


exp 


\-i ( r “ {Y{{ol + /3)/2),/3 - a)} = O (e < a - fd Y) . 

[ JAA] ) 


Using the rapid decay of with respect to 7 S F, we conclude that the spectrum of the 

matrix in (13.541) is at a Hausdorff distance of order e from the spectrum of 


exp{-i (He((Q: +/3)/2),/3 - a)} |B*| W eydp{i{9, a - (3)}/iim{0)d6 in Z^(r)^. (3.55) 

JE, 

The above matrix can be identified with a bounded operator in L^(T’)^ ^ [^^(r) 0 L‘^{E)]^ by 
taking the tensor product with the identity operator in L'^{E). Its N x N matrix-valued integral 
kernel is then given by: 

K'^{[x] + X, [x'] -I- x') := S{x — x') exp {—i (He(([x] -I- [x'])/2), [x'] — [x])} /I([x] — [x']). (3.56) 

This operator is isospectral with the matrix in (I3.55|) . hence its spectrum lies at a Hausdorff 
distance of order e from the spectrum of H^Ej{H'^). 

Next we compute the integral kernel of the operator Dp(/r'^) where /E{x,^) = — A^{x)), 

and we shall compare it with K^. The symbol is r*-periodic in the ^ variable. The N x N 

matrix valued integral kernel of its Weyl quantization acting on is given by: 

^^([x] + X, [x'] + x') := (27r)-‘' [ /r^((x + x')/2, (3-57) 

J X* 

Each ^ G X* can be uniquely written as ^ = [^] + 0 with [^] G and 9 G E^. We have 
/A{x,^) = /A{x,9) and using 1) = 1 together with the completeness relation in L^{E) 


2 ^ IK...1 


6{x — x') 
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we have: 


C^(x,x')= / d9{2Tr)-‘‘ V +^'0/2+ (N + M)/2,0) 

Je, -r»er, 

= 6ix - i')\E,\-^ [ - A,{x + ([:r] + [x'])/2)) 

J 

= 5{x - x')exp{-i {A^{x + ([x] + [x'])/2), [x'] - [x])}(3.58) 

Because E is bounded and A has all its first order derivatives bounded we obtain the estimate: 

(A,(x + ([x] + [x'])/2), [x'] - [x]) - (^.(([x] + [x'])/2), [x'] - [x]) = 0(e < [x'] - [x] >). 

Due to the localization properties of /I.y, we can replace A^[x + {[x\ + [x'])/2) with ^e(([x] + [x'])/2) 
in (13.5811 with the price of an error of order e in the norm topology. But then we obtain the kernel 
in (13.541) and the proof is over. 
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